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The purpose of this brief note is to present an infinite dimensional family 
of exact solutions of the incompressible three-dimensional Euler equations 

^tu + u ■ Vu + Vp = 0, V ■ u = 0. (1) 

The solutions we present have infinite kinetic energy and blow up in finite 
time. Blow up of other similar infinite energy solutions of Euler equations has 
been proved before The particular type of solution we will describe 

was proposed in (0). The Eulerian-Lagrangian approach we take (0) is not 
restricted to this particular case, but exact integration of the equations is. 
We will consider a two dimensional basic square Q of side L. The particular 
form of the solutions (0) is 

u(x, y, z, t) = {u{x, y, t), z-f{x, y, t)) (2) 

where the scalar valued function 7 is periodic in both spatial variables with 
period L and the two dimensional vector u{x,y,t) = {ui{x,y,t),U2{x,y,t)) 
is also periodic with the same period. The associated two dimensional curl 
is 

du2{x,y,t) dui{x,y,t) 
-(-'^'^) = ^^ 

This represents the vertical (third) component of the vorticity V x u of the 
Euler system and, using the ansatz (H) it follows from the familiar three 
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dimensional vorticity equation that the equation 



dtOJ + u ■ Voj = 70; (4) 

should be satisfied for the recipe to succeed. On the other hand, one can 
easily check that the vertical component of the velocity z^{x,y,t) solves the 
vertical component of the Euler equations if 7 solves the non-local Riccati 
equation 

|l + ^.V7 = -y + J(t) (5) 

where I{t) is a function that depends on time only. The divergence-free 
condition for u becomes 

V ■ M = -7. (6) 

Because of the spatial periodicity of u one must make sure that 

-f{x,t)dx = (7) 

Q 

holds throughout the evolution. This can be done provided the function I{t) 
is given by 

where 

\Q\ = [ dx = L\ 
JQ 

The velocity is determined from ui and 7 using a stream function ip{x,y,t) 
and a potential h{x, y, t) by 

u = V^ip + V/i (9) 

- A/i = 7, (10) 

-Atlj = uj (11) 

with periodic boundary conditions. 

The ansatz u{x,y,z,t) = {u{x,y,t), Z'y{x,y,t)) associates to solutions of 
the system (^, |, ||, ^, 0, |llj) in = 2 velocities u that obey the incompress- 
ible three dimensional Euler equations ([^, f^). The divergence condition (|^) 
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follows from |TOD. The compatibility condition Jqudx = is maintained 
throughout the evolution because of (|, ^. We consider initial data 

7(x,|/,0) = 7o(x,?/), uj{x,y,0) = uJo{x,y) (12) 

that are smooth and have mean zero Jq '-jodx = Jq uj^dx = 0. The solutions 
of the system above have local existence and the velocity is as smooth as 
the initial data are, as long as sup^ \'y{x,t)\dt is finite (the result can be 
proved (|^) following the idea of the proof of the well-known result (0))- 
We will consider the characteristics 

d X 

^ = uiX,t) (13) 

and denoting X{a, t) the characteristic that starts at t = from a, X{a, 0) = 
a, we note that, prior to blow up the map a i— >• X{a,t) is one-to-one and 
onto as a map from = R^/LZ^ to itself. The injectivity follows from the 
uniqueness of solutions of ordinary differential equations. The surjectivity 
can be proved by reversing time on characteristics, which can be done as 
long as the velocity is smooth. Our result is an explicit formula for 7 on 
characteristics 

where A{x,t) is the inverse of X{a,t) (the "back-to-labels" map) and the 
functions T(t), ^(t) and 0(r) are computed from the initial datum 79. More 
precisely we show 

Theorem 1 Consider the nonlocal conservative Riccati system (0, ^ 



[7q , \T]\ ). There exist smooth, mean zero initial data for which the solution 
becomes infinite in finite time. Both the maximum and the minimum values 
of the solution 7 diverge, to plus infinity and respectively to negative infinity 
at the blow up time. There is no initial datum for which only the minimum 
diverges. The general solution is given on characteristics in terms of the 
initial data 7(x, 0) = 7o(a;) by (see (^^) 

l{X{a,t),t) = a{T{t)) f .^^°Sf . . -mt))] 

\l + r(t)7o(a) / 
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where ^ 

'^^^ ^ -da\ I f —doS 



g (1 + r7o(a))2 J [Jq 1 + r7o(a) 



\Q\ Jq 1 + r7o(a) 
and 

^ = a(r), r(0) = 0. 
The function r(t) can 6e obtained from 

r/ie Jacobian J{a,t) = Det |^^^^| g'zi'en 6?/ 

1 fir da 1 ^ 



J(a,t) 



1 + r(t)7o(a) [ \Q\ Jq 1 + r(t)7o(a) 
The moments of 'j are given by 



p 



{^{x,t)fdx = {a{T)Y f , , -<P{rm Jia,t)da. 

Jq [1 + r(t)7o(a) J 

The blow up time t = T^ is given by 

Iff 1 , ( 7o(a) - mo \ ^ 



mo = min7o(a) < 0. 

Q 

We note that the curl uj plays a secondary role in this calculation and in 
the blow up. Indeed, the same formula and blow up occurs if ojq = 0, or if 
the curl u was smooth and computed in a different fashion than via (HI). 
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1 Solving on characteristics 

We will solve now the nonlocal Riccati equation on characteristics. We start 
with an auxiliary problem. Let (f) solve 

dr(i) + W ■ V</) = -0^ 

together with 

V • v(x, t) = —4>(x, t) + — — / 0(a;, T)dx. 

\Q\ JQ 

We will consider initial data that are smooth, periodic and have zero mean, 

/ 4>o(x)dx — 0. 
JQ 

We will also assume that the curl C — ^ — ^ obeys 

^ ax ay 

+ ^ • vc = (0 - 1^ <P{x, T)dx^ c- 

Passing to characteristics 

^-v{Y,r) (15) 

we integrate and obtain 

y(a,r),T) = 

1 + r0o(a) 

valid as long 

inf(l + T0o(a)) > 0. 
aeQ 

We need to compute 



\Q\ 

The Jacobian 

J(a, r) = Det 



dY' 
da 



obeys 

— = —h{a,T)J{a,T) 
dr 

where 

h{a, t) = (j){Y{a, t),t) - ^^{t). 
Initially the Jacobian equals to one, so 

J(a, r) = 

So 

J(a, r) = e/o^^(^)'^^ exp (- f log(l + s0o(a))rfs) 

JO as 

and thus 

J(a,r) = e/o^^^')'^' 



1 + r0o(a) 

The map a i— > l^(a, r) is one and onto. The change of variables formula gives 
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{x,T)dx = / (j){Y{a,T),t)J{a,T)da 

JQ 

and therefore 

\Q\ Jq (1 + r0o(a)) 
Consequently 

dr dr |Q| 1 + r0o(a) 

Because both sides at r = equal one, we have 



r) = e/;^W^^i-/^^^^4^da. (16) 



IQI Jq 1 + r0o(a 



-da (17) 



and, using (piq). 



^) = <! / 7^-^#7-^^4 ^ / . , . , da \ (18) 



(l + r0o(a.))^ J [JQl + T(f)o{a 
Note that the function 6{x, r) = 0(x, r) — 0(r) obeys 

(9^ 1 /• _ 

— + . V5 = + 2— / 52^3; _ 205, 
9r IQI 7q 
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We consider now the function 
and the velocity 

U{x,t) =e^^o^^'^'^%{x,T). 
Multiplying the equation of S by e^-^o ^(*)'^^ -^q obtain 



OT \Q\ J 



Note that 



V • [/ = -(7. 

Now we change the time scale. We define a new time t by the equation 



t(0) = 0, and new variables 



and 



dr 



-f{x,t) ^a{x,T) 



(19) 



u{x, t) = U{x, t). 
Now 7 solves the nonlocal conservative Riccati equation 



'di 



with periodic boundary conditions, 

u={-A)-^ [V^a; + V7 

and 



'dt 



+ u ■ Vu) — ^uj 



(20) 

(21) 
(22) 



The initial data are given by 

lo{x) = 5o{x) = 0o(a;). 
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Using (|T^ and integrating the equation ( |I9| ) we see that the time change is 
given by the formula 
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1 \ f f 1 , l + T(j)oia) 



'=m) JJo Ma) -Mb) TTT^""* 

Note that the characteristic system 

- = ..(A-.t) 

is solved by 

X{a,t) = Y{a,T) 
where Y solves the system ([T3|). This implies the formula 

7(X(a, t), t) = a(r) f t^^^ - 0(r)) (24) 

with 

a(r) = e2/o^<^W''^ (25) 



In view of (|17D, (|Tq) , (|23|) , we have obtained a complete description of the 
general solution in terms of the initial data. 



2 Blow up 

Consider an initial smooth function 70(a) = 0o(c^) and assume that it has 
mean zero and that its minimum is tuq < 0. As it is evident from the explicit 
formula the blow up time for (f){Y{a, r), r) is 

1 

mo 

and (j){Y{a,T),T) diverges to negative infinity for some a, and not at all for 
others. This of course does not necessarily mean that 7 blows up in the same 
fashion. Let us discuss the simplest case, in which the minimum is attained 
at a finite number of locations ao, and near these locations, the function 0o 
has non-vanishing second derivatives, so that locally 

0o(a) > mo + C\a - aop 
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for < |a — ao| < r. Then it follows that the integral 

If da 

\Q\ 7 e2 + 0o(a) - mo 

behaves like 



\Q\ J + (j)o{a) - mo ] V V e 

for small e. Taking 

r 

we deduce that, for these kinds of initial data 



1+^ 



r* — r 



For the same kind of functions and small (r^, — r) , the integral 

Ma) . C 



\Q\ Jq (1 + TMa)y T.-T 

and t(r) has a finite limit t —>■ as t —>■ t^:. The average 0(r) diverges to 
negative infinity, 



ir)' ^ 



n — T 



log 



C 



— r 



The prefactor a becomes vanishingly small 

a(r) ~ (log(r* - r))"^ 

and (^) becomes 

0o(a) 



7(X(a,t),t) ~ (log(r,-r)) 



-2 



1 + r</)n(a) 



If the label is chosen so that 0o(a) > then the first term in the brackets 
does not blow up and 7 diverges to plus infinity. If the label is chosen at the 



9 



minimum, or nearby, then the first term in the brackets dominates and the 
blow up is to negative infinity, as expected from the ODE. From the equation 

m 

(a(r)) -^dr = dt 

it follows that 

- t ~ (r, - r) (l + log 



and the asymptotic behavior of the blow up in t follows from the one in r. 
We end by addressing a question that was at some point raised by numerical 
simulations: can there be a one-sided blow up? From the representation (|^) 
of the solution it follows that 

M{t) > -0(r)e'/o^^W'^* 

holds for M{t) = sup^7(x,t). If one would assume that, up to the putative 
blow up 

M{t) < C 

with some fixed constant C, then it would follow that 

dr 

and integrating between r and that 

This in turn would imply that T^, = oo and therefore no blow up for 7 can 
occur in finite t. So the answer is that for no initial datum can there exist a 
one sided blow up for 7. 
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